General solution for Hamiltonians with extended cubic 

and quartic potentials * 

C. Verhoeven, M. Musette and R. Conte 
10 January 2003 



Keywords. Henon-Heiles, Hamiltonian systems, separation of variables, nonlinear 
equations, hyperelliptic functions, soliton equations. 

Abstract 

We integrate with hyperelliptic functions a two-particle Hamiltonian with 
quartic potential and additionnal linear and nonpolynomial terms in the Liou- 
ville integrable cases 1:6:1 and 1:6:8. 

1 Introduction 

The generalized Henon-Heiles Hamiltonian 

H = hpl^ + + ciX^ + C2Y^) + aXY"^ - ^X^ + fiY-\ /x arbitrary (1) 

is known to be Liouville integrable |B 121 EI for three sets of values of {b/a, ci, C2) 
and to be related [5 to the stationary reduction of the following fifth order soliton 
equations: KdVs, Sawada-Kotera (SK) and Kaup-Kupershmidt (KK). 
The canonical transformation [S] between the SK and KK cases for /i 7^ allowed us 
[Sj to define the separating variables of the Hamilton- Jacobi equations and to derive 
the general solution of the equations of motion with hyperelliptic functions. 
The two-particle Hamiltonian with quartic potential 

i7 = 1(P| + _ 1 (qJ^2 ^ ^y2) ^ ^^4 ^ BX'^Y'^ + AY^, (2) 
is known to be Liouville integrable in four cases [3 |H1 [U] : 
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4. A : B : C = 1 : 12 : 16,a = 46. 

The extension of Q to include some linear and non polynomial terms which preserve 
the Liouville integrability jlUl [U] and confirm their connection with some integrable 
soliton equations has been considered by various authors |1H I12j . 
The case 1 with extra inverse square terms is equivalent to the travelling wave 
reduction of the Manakov system of two coupled NLS equations. It corresponds 
to the particular case for two particles of the extended Garnier system ^1] which 
has been integrated by Wojciekowski The case 4 can be associated jH] with a 
coupled KdV system possessing a fifth order Lax pair . 

The authors of |17| |^ have showed how cases 2 and 3 with some extra linear and 
inverse square terms are equivalent to the stationary reduction ^ = x — ct of the 
Hirota-Satsuma system |18j : 

_ 1 

ft — 2^^^^ ~^ "^f^x GsSx, St — Sxxx 3rSj; (3) 
and an other coupled KdV system 

ft = -Wfxxx + iff XX + 3/2 - 3/2/x + Qf9x + 6gfx)^ 

gt= ji2gxxx + i'^ggx + Qfgxx + i'2gfxx + i8fxgx-Qf'^gx (4) 

'i'ifxXXX + if f XXX + ^8fxfxx ~ 6/ fxx ~ ^ffx)'' 

In the present paper, using the canonical transformation 5^ between these two ex- 
tended cases, we show how the method we followed for the SK and KK cases of Q 
to define the separating variables of the Hamilton-Jacobi equation can be applied 
here. Then the equations of motion associated with cases 2 and 3 of © with extra 
linear and inverse square terms are also integrated with hyperelliptic functions. 
In section 121 we consider the SK and KK integrable cases of the Henon-Heiles Hamil- 
tonian for /u 7^ to explain the method we applied in JB^ for integrating their 
equations of motion. In section |31 we consider the Hamiltonian Q with extended 
quartic potential in case 2 and 3 and recall the canonical transformation 5_, which 
relates the case 1 : 6 : 1 of ((21) (with extra terms of the form aX~'^ + (3Y~'^) and the 
case 1:6:8 (with extra terms like ^yX + 5Y~'^). 

To the best of our knowledge, the results of sections and 13 are new. In section HI 
we explicitly integrate the equations of motion of the case 1:6:8 for 5 7^ 0, 7 = 0. 
In section ISl we integrate the 1:6:1 case with a = (3 ^ and we use the canonical 
transformation to transport those solutions to the 1:6:8 case for (5 = 0, 7 7^ 0. 

2 Integration of the cubic Henon-Heiles Hamiltonian 

Let us consider the SK and KK Henon-Heiles Hamiltonians with their second inte- 
gral, denoted as K2 and /cl , and their equations of motion ^ : 

SK: - = -1, ^ci=C2, « = -^! U = X + C2, V = Y, c = C1C2 (5) 
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Kl = Kl, + ^f,U + fi^, (7) 

i^2,o = -2PuPv - f/V - — + cF, (8) 

[/" = -i(F^ + [/2)-|, V" = -UV + ^^, (9) 

KK: - = -16, ci = 16c2, a = \, c = ciC2, ii = X + 2c2, V = HlO) 
a 4 



^ = A:i = ^(p'+p') + 7^2 + 1^3-0^+^4' (11) 
2 4 J 2 



:,2 _ 7^2 I M„, , o,.!'^ I 



^2 = ^^2,0 + ^^+2^^ + ^' (12) 

^2,0 = - 7^^^ - ^2^^^"^ ^^^^^ ~ gP^P^'''^^ + Y^^^ (13) 

n" = -lt;2-47x2 + c, ^" = -in^; + 4, (14) 
4 2 ^3 

which are equivalent to the stationary reduction = x — ct of the two integrable 
PDE's: 

SK: Ut + {a,,,.,, + 5UU,, + ^U% = 0, (15) 

20 o 9 

KK: ut + (u^^^^ + lOnn^.^. + y + 30<)^ = 0. (16) 



Both equations are connected to another fifth order integrable PDE 

wt + (w4x - 5wxWxx - ^w^Wxx - 5wwl + W^)x = 0. (17) 
by the Miura transformation 

U = Wx — w"^, u = —Wx — 2^^' (1^) 

Solving the stationary reduction of (|17|) for w in function of V, V' and v, v' and 
defining the following expressions: 

A' = -/i, r = 6{VK2,o + XPu), 

+3A\v^ + ISX'^v'^pI + UX^vp^ + 3A^), (19) 
one obtains the following canonical transformation [Sj l2Uj : 



Pu = ^{^Pv + ^UV^Pv - PuV^) 
1 / A\ /- . y 



Pv 



2Pv + -]^-\^, (21) 



4fH Vr 
- + — -u, ^ = Tn 



U = -^- + -2 -u, V = ^^, (22) 



Pu = \{l2pl + Guv'^py - v^pu) 



3A 

+^i2uv^ + 12v^pI + 12Aup^ + 4A2), (23) 
v" 



= -l^{p, + ^^^^ + X—. (24) 



For = 0, we introduce the expressions (|22j) - (|23)) for U,V, Pu, Py in the variables 
which separate the SK Hamiltonian 



Q,=U + V, Pi = \{Pu+Pv)^ 
Q2 = U-V, P2 = \{Pu-Pv)- 



(25) 



This defines the change of variables |2 

- k2fl 1 

— u, p^ = — 

2v^ 



qi = -Q— — 5^ -u, pi = -^{12pl + Guv'^py - v^pu - 12p^k2fi) 



q2 = -6 ^" ~^2^'° -u, p2 = -^{I2pl + 6uv^py - v^pu + 12p„A:2,o)' 



(26) 



that we apply on the KK Hamiltonian taking account that for fi ^ 0, A;2,o is 
no more a constant of motion. Therefore, one has 

H^h=pl+pl + ^^{qf + 4) -l{qi + 92) + (27) 

k2fi = 2{pI - pI) + \{ql- ql) - ^(^2 - qi). (28) 
2 

/ „ , ^J' {<ii-q2)pi ,r,n\ 

/ o f^iqi- 92)P2 .on\ 

92 = 2P2 - p . (30) 

U «,2,0 

In this new setting of coordinates, it is important to note that, introducing the 
expression 

fiqi,Pi) = 2pi + \ql - for i = 1,2, (31) 
2 
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in (|27|) and (|28j) . the Hamilton- Jacobi equation is separated for fx arbitrary. 

In the variables (|2H|) . the second invariant fc^ can be written in two equivalent ways: 

or ^2 = -392+(^A:2,o-Y^^J , (33) 

which, after the elimination of fi{qi — (?2)/fc2,o between (|27|) -(|32 |) and (|27]) ~ (j33|l . 

become 

fci = -|gi+(^-4p?-| + cgi + 2fci^ , (34) 

or kl = -i^q2+(^4pl + ^-cq2-2h^ ■ (35) 

Then we can eliminate pi between (|3l|) and (|29|) . and p2 between (|35|) and (|30|l to 
obtain 



l[ = \l2k,-i + cq, ,lkl + g,, I 1 + ^ _ I , (36) 



'/^ - J2fci - I + + + ^l2\l- ^ , , '^^ ''^^ I • (37) 



fci + §92 + ./fci + §gi 



For n ^ 0, setting: 



si = y3^ + (/i, s2 = -y3^ + (72, (38) 



and defining 



P(.) =2^1-1 - ) ' + c - sj) - ^s, (39) 



the equations (|5^ -H37 |) become: 



= yZM, 4 = _ (40) 

Si - S2 Si - S2 



They can be solved by inverting the hyperelliptic integrals 



''I ds r2 ds , 

+f^=f + ,„ (42) 



00 



VPis) Joo VPi^ 
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and the general solution for the KK system is 

2 



"= 2 



\(4 + 4) + ->^i-l('-^) . v' = ^. (43) 

2 /X 2 VSI + 82/ Sl + S2 



For ^ = from ()36() - (|37p and ()26() . one easily recovers the known solution |21j 
expressed with Weierstrass elliptic functions: 

3 / pi - P2 \ ^ 2 _ ^2,0 



2 V Pi - P2 / P2 - Pi 



where p, satisfies the equation 



pf = 4p3_52P^-5r. ^ = 1'2 (45) 
52 = ^, = - — (2/21,0 -/c2,o), 9^ = -^44 (2^1.0 + ^2,o)- 

In the SK case, from ()43() and the canonical transformation H22() . the general solution 
writes 

u = V^{s[ + 4) + + siS2 + si - 

^2 = -2^{si + S2)(S1S1 + S2S'2) + 2{si + ^2)^ (sf + - ^) , ^^^^ 

which, in the particular case /i = 0, is merely 

C/ = -6(pi + p2), F = -6(pi-p2)- (47) 

3 Canonical transformation between extended 1:6:1 
and 1:6:8 Hamiltonians 

Let us consider the quartic Hamiltonian @ with extra linear and inverse quadratic 
terms in the two following cases jH]: 

1:6:1 H = Ki ^ ^{Pi + P}) - j^{U' + 6f/V= + V') 

+ + + (48) 

Kl = Kl, - PU' - aV' + 4(a^ + (3^) + - 4c(a + /?), (49) 

K2,o = 2PuPv - \UV{U'^ + V'^ + 8c), (50) 

U' = Pu, V = Pv^ (51) 
U" = + luv' + cU + ^ V" = + lu'V + + (52) 



6 



1:6:8 H = h = ^ipl+pl) - ^{8u^ + Gu^v^ + v^) 
2 lb 



|(4n2 + t;2)-^n + ^, (53) 



+A-f{vpuPv - upl) - 475^ - 26c - Ac-yv'^u, 
It' = pu, v' = (56) 



(54) 
(55) 



3 1 3 (5 

u" = 2u^ + -uv'^ + 4cu + -y, v" = -v^ + -u^v + cv + (57) 

4 4 4 

where i^i, and ki, ^2 are the constants of motion of both systems. 
The equations of motion are respectively equivalent to the stationary reduction 
^ = X — ci of the PDE's © and Q with the correspondence: 

r = -hu^ + V^ + 4c), s = l{U^-V^), (58) 
4 o 

f = u, g = -^iv^ + 2u^ + 2p^ + 4c)- (59) 

Both PDE's © and Q possess a fourth order Lax pair with scattering operators L 
and L which factorize as follows |5l I17j. 

1:6:1 L = {d^ + r + s){dl + r - s) 

= {dx-vi){dx + vi){dx + V2){dx-V2)^ 

1:6:8 Z = {d^ + fd, + f, + g){d, - fd, + g) ^^^^ 

= {dx + Vl){dx + V2){dx - V2){dx - Vi), 

and yield the Miura maps: 

1:6:1 r = ^{vi^x - V2,x - vf - V2), s = ^{vi^x + V2,x - vf + V2)^ (62) 
1 : 6 : 8 / = -ui + i;2, g = V1V2 - vi^x^ (63) 
while vi and V2 are solutions of the system of PDE's 

Vl,t = l(-2ui,a;x - GV2,xx " 12^2^2,2: " 12flW2,x " l2viV^ + 4:vf)x^ ^ ^ 

(64) 

V2,t = |(-2W2,xa; " Qvi^xx + 12flfl,x + 12f2'Vl,x - 12i;2Wi + 4:V^)x- 

From the relations (|62|) - l|63j) . the transformations (|58|) -(|59 |) and the equations of 
motion ((^ - H52|) . ((^ - H57|) . the stationary reduction of 1)64^ can be solved for vi 
and V2 in function of the canonical variables of the two extended systems: 

1:6:1 vi = -— + -^, V2 = — + -^, (65) 



u Pv V —6 u py \/—d 
1:6:8 vi = - -\ \ 5-, f 2 = tt 5-' (66) 
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such that, defining 

a = -K\, P = -kI, (5 = -(k2 - Ki)^, 7 = i(Ki+K2)i (67) 

^ = -\iU' + V') + ^{PuPv + -2^^ - - ^) - 2- (68) 
= -4( T 2p„ + + - 44 ± 4^ + 8(^2 - ki)^ 

^4(^2-^1)4-4^^^^^ + 4c), (69) 

the canonical transformation between the 1:6:1 and 1:6:8 cases becomes jHj 

Pu = -^{V -[/) + _- — - 2^1^ - 2k2^ + [74 - 174' (70) 
..270 ,„.^(_^_^.-_-)_--. 

J/^yrr, P,.4vrr(|.^-^i^KWr-. (72) 



^ = 2^' ^^ = 2^U"v + ^;^j-^"^^- ^^^^ 

4 General solution of the 1:6:8 Hamiltonian for 7 = 

For a = P = 0, one introduces the transformation (|72|) - (|7:'{j) in the variables 

Qi = i(c/ + y)2, p, 



IPu + Pv 



2 U + V (jA\ 
-,yU-V) , P2-2-^73T^' 



which separated the 1:6:1 Hamiltonian 1)48^ in the polynomial case. This defines 
the following change of variables: 



Pi 

Q2 
P2 



4^^'°+^' 4c, 
2 




(75) 




>k2,o) + 8cui;^ 


(76) 


16v{vpu - 2upv) 




2 




(77) 


—SvpvPu + 8upl + 2f + u{v^ — 5 


)fc2,o) + 8cuu^ 


(78) 


16-u(vp„ - 2upv) 
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that we apply on the Hamiltonian system (|53|) for K1K2 7^ 0, i.e. when A;2,o is no 
more a constant of motion. This transformation restricted to K2 = —Ki {S = — 4/^1 , 
7 = 0) defines the coordinates which separate |lllll2j the Hamilton- Jacobi equation. 
Indeed, one has 

H = h = 3^(32gip? -ql + 32q2pl - ql - Sc{qi + ^72)) - (^9) 



with 

^2,0 



^ (32gip2 _ ql - 32q2pl + ql + 8c(g2 - qi)) ' (80) 



91 = 4gipi + 2/ti -2 ) (81) 

92 = ^q2V2 - 2/€i -2 > (82) 

'*^2,0 



such that defining 

f{qi,Pi) = ^'^qiPi - qf - 8cqi, i = i, 2 

we have the following separated Hamilton-Jacobi equation: 
f{qi,Pif-f{q2,P2f-^2Kl{qi-q2) = I6ki {f{qi,pi) - f{q2,P2)) , Pi = (83) 



dS 



dqi 

In analogy with the Kaup-Kupershmidt case, we can write the second invariant in 
two equivalent ways: 

kl = -2nlq^ + ( k2,o + (84) 



2k2,o , 

or kl = -2Klq2 + (k2,o - ^%^) (85) 

which allow us to ehminate Acf (gi - Q'2)/^2,o between ((TU-dHH) and Next, 
we eliminate pi and p2 between those two resulting expressions and the equations 
(jHT|) . (jS21) and obtain: 



1[ - J'-k + 4c + 4^ + ^J^~^A I + ^ (86) 



2 qi qi 



(^^kj + 2K,lqi + ^kl + 2^292 



. J? + 4c + 4^ ^ 1 - (87) 



2 g2 92 



v/Ai + 2«^ + yAi + 2fi^ 



For Ki = the general solution of equations (|86]) - (|87|) can be expressed in terms of 
Weierstrass elliptic functions: 

8 . 4 2 fei A;2 c , , , 32 2^\ 

91 + 3C = 8p(e - 6, - Y - + /C2 - —c^)) = 8pi' 

8 / ^ 4 n /ci A;2 c , , , 32 9, \ , , 

92 + 3C = 8p(e - 6, - y + f , Y^C^h -k2- -c')) ^ 8p2- (89) 
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For Ki^O setting in (|86|l -(l871): 



si = W^^2 + 9i and ^2 = -W^fc2 + 92, (90) 



and defining 



we obtain the equations, 



Si - S2 Si - S2 

which are solved with the inversion of the hyperehiptic integrals 



''I ds ds , 

= + / = ^3, (93) 



oo 



Pis) -"^ \ Pis 



''I sds ^2 sds ^ , , 

= + / = e + ^4- (94) 



oo 



Pis) -"^ \ Pis 



Therefore, for 6 arbitrary, 7 = 0, the general solution of the 1:6:8 Hamiltonian is 
defined with symmetric combinations of si, S2 



w = — (si + So H 0+ 4c, V = ) (95) 

^ ^' VS1 + S2/ S1 + S2 ' S1+S2 ^ ' 



and is a single- valued function of the complex variable ^. For ki = it degenerates 
into 

u = -4(pi + P2) + — + ttt; r - :tC, v = — • (96) 

5 General solution of the extended 1:6:1 Hamiltonian 

Now, we again use the canonical transformation H72p - H73() to find the general solution 
of the 1:6:1 Hamiltonian. 

For Ki = K2 = 0, we introduce the expressions in ((7 ^ - H73p and obtain the 
solution of the equations of motion 

U^ + V^ = 8(pi + p2 - \c), UV = 4(pi - p2) • (97) 
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Next, for af3 7^ 0, with the restriction a = /?, we obtain the expressions 

U^ + V^ = -2V2{s[ + 4) + 2{sl + sl + S1S2) - (98) 
U^V^ = -2\/2(si + S2){sis[ + - 2ki) 

+2(si + S2)'(^5? + si-|^ + 4c^, (99) 

which depend on symmetric combinations of si,S2 and therefore are single- valued 
functions of S,. 

Since the Hamiltonian (|48() for a = /3 is an even function of ki, we apply on H98() - (|99() 
the canonical transformation (|7()|) - (|71j) for ki = K2, and obtain the general solution 
of the 1:6:8 Hamiltonian in the extended case 5 = 0. 

6 Conclusion 

Romeiras described a procedure which relates the extended cubic potential 
V3{u,v) = V3{u,v) + fiv-'^/2 in the KK case with the extended quartic potential 
t4(n, v) = V4(u, v) + 5v-^/2 in the 1 : 6 : 8 case. This explains why the combination 
u'^ + v'^/2 + 4c of the solutions is identical to the solution H43() . The question 
remains open to extrapolate the present results to the case a 7^ /?, i.e. 7(5 7^ 0. 
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